Novel Quantum Criticality in CeRu 2 Si 2 near Absolute Zero Observed by Thermal Expansion and 

Magnetostriction 
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We report linear thermal expansion and magnetostriction measurements for CeRu2Si2 in magnetic fields up to 
52.6 mT and at temperatures down to 1 mK. At high temperatures, this compound showed Landau-Fermi-liquid 
behavior: The linear thermal expansion coefficient and the magnetostriction coefficient were proportional to 
the temperature and magnetic field, respectively. In contrast, a pronounced non-Fermi-liquid effect was found 
below 50 mK. The negative contribution of thermal expansion and magnetostriction suggests the existence of 
an additional quantum critical point. 

PACS numbers: 71.10.Iff, 71.27.+a, 75.80.+q 



A quantum phase transition (QPT), a continuous order- 
disorder transition at zero temperature, is driven by quan- 
tum fluctuations and differs essentially from a classical phase 
transition, which occurs due to thermal fluctuations of an or- 
der parameter A 4/-based heavy-fermion (HF) system is 
suitable for studying quantum critical behavior and has at- 
tracted experimental and theoretical attention. The ground 
state in this system, determined by the competition between 
the on-site Kondo interaction and the inter-site Ruderman- 
Kittel-Kasuya-Yosida interaction, can be tuned by physical 
parameters such as pressure, elemental substitution, and mag- 
netic field. Furthermore, near a quantum critical point (QCP), 
quantum fluctuations cause a pronounced deviation from the 
Landau-Fermi-liquid (LFL) behavior, leading to what is re- 
ferred to as non-Fermi-liquid (NFL) behavior. In the tradi- 
tional picture J2, S 01, quasiparticles retain their itinerant 
character and form a spin-density wave (SDW) type of an- 
tiferromagnetic (AFM) order. This itinerant model has suc- 
cessfully explained experimental critical phenomena close to 
QCPs [5]. However, recent experiments have shown that this 
model fails in some HF systems (a, [3]. Thus, a new type of 
QCP associated with the decomposition of quasiparticles has 
been proposed lalSl]. 

It has recently been shown that the thermal expansion co- 
efficient a is more singular than the specific heat C in the 
approach to a QCP. Accordingly, the Griineisen parameter 
r oc a/C diverges as temperature goes to zero near the QCP 
and the associated critical exponent of T can be used to deter- 
mine the type of QCP 11011 . In addition, the Griineisen param- 
eter and the thermal expansion coefficient exhibit a character- 
istic sign change close to the QCP [11]. Thus, the Griineisen 
parameter and the thermal expansion coefficient are widely 
used as excellent tools for the investigation of QCPs in vari- 
ous systems, for example, HFs Il2lll3ll . spin-ladders 11411 . and 
perovskite ruthenate [15]. 

The intermetallic compound CeRu2Si2, which has a 
ThCr2Si2-type crystal structure, is well known as a typical 
HF compound, and exhibits the LFL state with an electronic 
specific heat coefficient 7 = 350 mJ/K 2 mol below the Kondo 



temperature Tk = 14 K |3, 171. Neutron-scattering mea- 
surements represent short-range time-fluctuating AFM corre- 
lations described by three incommensurate wave vectors even 
below Tk flUl . The elemental substitutions Cei_ K La a; Ru2Si2 
and Ce(Rui_ a; Rh 2 ;)2Si2 induce an incommensurate SDW 
ground state at a finite transition temperature. With decreas- 
ing x, the transition temperature decreases to zero at the fi- 
nite critical concentrations x c = 0.08 in Cei_ a; La :E Ru2Si2 
and x c = 0.03 in Ce(Rui_ ;c Rh 2; )2Si2 dlM]. Since these 
substitutions induce an expansion of the lattice parameter, 
CeRu2Si2 is regarded as an example of a pres sure-driven QPT 
with a critical pressure p c = —0.3 GPa [21]. In the elemen- 
tal substituted system close to the critical concentration, sev- 
eral critical phenomena down to 0.1 K have been explained 
by the 3D spin fluctuations of the itinerant model 
In contrast to this traditional picture, our previous experi- 
mental results for undoped CeRu2Si2 at small magnetic fields 
and ultralow temperatures have shown the NFL behavior, that 
is, an ac susceptibility peak, and a temperature-independent 
but field-dependent ac susceptibility and magnetization be- 
low the peak temperature of ac susceptibility, Tp B23I1 . The 
critical exponent of the susceptibility below 10 mK could not 
be explained by the itinerant model, and the suppression of 
susceptibility by the magnetic field showed a resemblance 
to the magnetic-field-tuned QCP observed in YbRfi2Si2 and 
CeQii-^Aua; In this Letter, we report thermal ex- 

pansion and magnetostriction measurements for CeRu2Si2 in 
magnetic fields up to 52.6 mT at temperatures down to 1 mK. 
We present the material's unique critical properties and dis- 
cuss the existence of an additional QCP. 

A single crystal of CeRu2Si2, grown along the c axis by 
the Czochralski pulling method using the starting materials Ce 
(99.99% pure), Ru (99.99% pure), and Si (99.999% pure), was 
cut into a cylinder, 3 mm in diameter and 5 mm in length. The 
sample was cooled with a copper nuclear demagnetization re- 
frigerator and a 3 He- 4 He dilution refrigerator. The thermal ex- 
pansion and magnetostriction along the c axis were measured 
by the capacitance method, with applying the magnetic field 
in the same direction. The bottom of the sample was glued to 
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the sample cell made of oxygen-free high-conductivity copper 
using Arzerite VL-10 silver paste for thermal conduction, and 
the top was epoxied to the 5-mm-diameter copper capacitance 
plate using Stycast 2850GT for electrical insulation. The dis- 
tance between this plate and another fixed capacitance plate 
was set to be 18 /im. The change in length AL of the sample 
was derived from the capacitance change. The sensitivity of 
AL/L was 10 -11 by the home-made capacitance bridge that 
was used. The thermal relaxation time observed in the exper- 
iment was less than several minutes around 10 mK. The mea- 
surements were done by continuously changing temperature 
and magnetic field at a constant rate from 0.4 to 4 mK/hour 
and 5.6 mT/hour, respectively. We observed no hysteresis 
in temperature and field sweep directions. Since the thermal 
expansion and magnetostriction of CeRu2Si2 are 4 orders of 
magnitude larger than those of Cu below 10 K 1124 12511 . the 
change in length of the copper holder was negligible. 

Figure 1 shows the temperature dependence of the lin- 
ear thermal expansion AL(T)/L in magnetic fields up to 
52.6 ml As shown in the inset of Fig. 1, the linear thermal 
expansion shows T 2 dependence, varying as AL(T)/L = 
aT 2 /2 in all measured fields above 50 mK. This tempera- 
ture dependence indicates that the system is in the LFL state, 
and the coefficient a = 5.60 x 10~ 6 K~ 2 obtained agrees 
with the results of Lacerda et al. [26]. However, with de- 
creasing temperature, AL(T)/L deviates from the T 2 de- 
pendence in all measured fields. The deviation temperature 
T* is about 50 mK at 0.016 mT and is suppressed slightly 
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FIG. 1: (color online). Linear thermal expansion AL(T)/L 
of CeRu2Si2 along the c axis at different applied fields, plot- 
ted as a function of T 2 . The data are subsequently shifted by 
lxlO~ 9 .The solid lines represent AL(T)/L = aT 2 /2 with a = 
5.60 x 10 -6 K _2 .The inset shows the temperature dependence of 
AL(T) / L at high temperatures. 



with increasing magnetic field. In a similar temperature 
range, the susceptibility shows a Curie-like increase from 
the temperature-independent Pauli paramagnetic susceptibil- 
ity 12311 . Thus, we believe that T* is the characteristic tem- 
perature of crossover from the LFL to the NFL state, and 
AL cr (T)/L = AL(T)/L - aT 2 /2, shown in the inset of 
Fig. 2, is the critical contribution accompanied by the NFL 
behavior. The critical contribution of the linear thermal co- 
efficient a a {T) is obtained by the temperature derivative of 
AL CI (T) I L. The temperature dependence of a cl -/T, shown 
in the main frame of Fig. 2, differs from 1 / \/T dependence, 
which is expected from the 3D spin fluctuations based on the 
itinerant model |_10D - In particular, the sign of a cr remains 
negative down to the lowest temperature, in contrast to the 
ordinary pressure-driven QPT of HF systems. 

In the control parameter vs temperature phase diagram with 
a symmetry-broken phase extended to a finite temperature 
from a QCP, the isentropes have a minima on the phase bound- 
ary and at the QCP [11]. In the case of a pressure-driven QPT, 
the sign of a cr is the same as that of the slope of the isentrope 
and changes on the phase boundary, according to the Maxwell 
relation (dV/dT) p = - {dS/dp) T . Most Ce-based NFL ma- 
terials have an AFM transition line at the low pressure side 
of the QCP in the p-T phase diagram. Thus, a a is positive 
when the pressure is higher than the phase boundary and neg- 
ative only inside the AFM ordered state. In the p-T phase 
diagram of CeRu2Si2, the AFM transition line ends the AFM 
QCP at -0.3 GPa from the substituted results lEltl . and no 
phase transition has been observed at ambient pressures down 
to 170 /iK M23I1 . Thus, undopedCeRu2Si2 at ambient pressure 
apparently stays on the high pressure side of the AFM QCP, 
and a cr should be positive down to the lowest temperature. 
One possible explanation for our negative critical contribution 
a CI is that an additional QCP exists in the region of higher- 
than-ambient pressure. Assuming that the undoped CeRu 2 Si2 
at ambient pressure is located closer to this additional QCP 
than to the AFM QCP, the isentropes have a negative slope 
owing to the isentrope minima of the additional QCP, and thus, 




20 30 40 
7~(mK) 



60 



FIG. 2: (color online). Critical contribution of the linear thermal ex- 
pansion coefficient a a (T) of CeRu2Si2, plotted as a a /T vs T. The 
inset shows the critical contribution of the linear thermal expansion 
AL CI (T) / L as a function of temperature. 
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a cr is negative down to the lowest temperature. 

Magnetostriction is a good measure of QCPs, since the 
magnetostriction coefficient A is related to the pressure 
derivative of magnetization via A = d[AV(B) /V]/dB = 
— (dM/dp)T,H on the basis of the Maxwell relation, where 
AV(B) /V is the volume magnetostriction and M the magne- 
tization per unit volume. Figure 3 shows the isothermal linear 
magnetostriction coefficient A(B) = L~ 1 d[L(B)—L(0)]/dB 
as a function of the magnetic field. Above 50 mK, A varies 
linearly with magnetic field, having a linear field coefficient 

A = dk/dB 9.6xl0~ 7 T~ 2 . In magnetic fields higher 

than those used in our measurements, A has shown a linear 
field dependence with a positive coefficient below 6 T and de- 
creased to roughly zero as the field decreased to zero 026112711 . 
The linear field coefficient can be estimated to be A ~ l.Ox 
10 -5 T -2 in the field range between 2 and 6 T at 70 mK 
112811 . Thus, the linear field coefficient is expected to change 
from negative to positive in a higher field. The negative sign 
with small absolute value of A observed in our low-field mea- 
surements implies that magnetization increases with pressure, 
which is expected given the quasiparticle mass enhancement 
in the approach to the QCP. Thus, the negative magnetostric- 
tion is consistent with the existence of an additional QCP at 
a higher-than-ambient pressure, as suggested by the negative 
thermal expansion coefficient. 

The additional QCP indicated by the thermal expansion and 
the magnetostriction at millike lvin temperatures in CeRu2Si2 
has a different cause than the AFM QCP. CeCu@ also shows a 
novel behavior at ultralow temperatures. Because the AFM 
transition temperature is suppressed to zero at the critical 
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FIG. 3: (color online). Isothermal linear magnetostriction coefficient 
A(B) = L _1 d[L(B) - L(0)]/dB of CeRu 2 Si 2 as a function of 
magnetic field. The solid lines represent the linear field dependence. 
For clarity, the data are shifted down by 1 x 10 T _1 at each temper- 
ature increment. Arrows show the crossover field to the linear field 
dependence. 



concentrations x c =0.1 in C&Cm§- x Axi x and x c = 0.09 in 
CeCug^Ag^ HI [ME 1], CeCu6 at ambient pressure is lo- 
cated on the high pressure side of the AFM QCP. In con- 
trast, CeCu 6 exhibits an AFM order at 2 mK, accompanied by 
the susceptibility peak and the jump in the thermal expansion 
coefficient, following the increase of susceptibility from the 
Pauli paramagnetic susceptibility below 100 mK 13211 . This 
behavior has been explained by AFM spin fluctuations. In 
the case of CeRu2Si2, the observed NFL behavior in suscep- 
tibility, magnetization, and a a - below 50 mK can hardly be 
explained by spin fluctuations based on the itinerant model, 
and the negative sign of a cr and A leads to the existence of an 
additional QCP. 

We will discuss the possible origin of the additional QCP. 
Yamagami and Hasegawa have proposed five Fermi surfaces 
on the basis of the band structure calculation 13311 . Indeed, 
the following four Fermi surfaces have been observed ex- 
perimentally: a large hole surface with 200 mo, a multicon- 
nected electron surface that has crossed arms with 10-20 mo, 
and two small ellipsoidal hole surfaces with effective mass 
1.5 mo and 1.6 mo, respectiverly, where mo is the free elec- 
tron mass 1134 13511 . The well known HF behavior around 1 K 
is mainly ascribed to the hybridization between the conduc- 
tion electrons and 4/ electrons (cf hybridization) on the large 
hole Fermi surface because of its large effective mass. We 
speculate that the effect of the cf hybridization between 4/ 
electrons and other Fermi surfaces with a small effective mass 
("small" Fermi surface) is not strong enough to contribute to 
the quasiparticle mass enhancement around 1 K, but this resid- 
ual hybridization becomes dominant with decreasing temper- 
ature. Owing to the spin fluctuations, the cf hybridization 
at the small Fermi surface further enhances the quasiparticle 
mass and the NFL behavior and may induce the additional 
QCP. 

The nature of the ground state in a magnetic field has been 
observed in the magnetostriction at lowest temperatures. As 
can be seen in Fig. 3, below 10 mK, a linear field dependence 
of A, indicative of the LFL state, is observed above 30 mT. 
However, A deviates from the linear field dependence below 
30 mT, which suggests the NFL state. The crossover field 
between the NFL and LFL states, indicated by the arrow in 
Fig. 3, decreases as the temperature is reduced. Such field- 
tuned crossovers are observed in many HFs. In the case of 
CeRu2Si2, the spin fluctuations due to the residual hybridiza- 
tion with the "small" Fermi surface that dominated below 
50 mK cause A to deviate from the linear field dependence in 
the low-field region, which agrees with the magnetization en- 
hancement observed in similar magnetic field and temperature 
regions [23]] ■ In the fields above the crossover field, the mag- 
netic field suppresses the spin fluctuations, which decreases 
the slope of A. 

Finally, we address the critical contribution of the 
Griineisen parameter, which is an important physical quan- 
tity characterizing quantum critical phenomena, as discussed 
in the introduction. The critical contribution of the Griineisen 
parameter can be expressed as r cr = V m a^ r / nC a , where 
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V m is the molar volume, the critical contribution of the 
volume expansion coefficient, k the isothermal compressibil- 
ity, C cr the critical contribution of the specific heat. Since 
the quadrupolar split of Ru nuclei prevents observation of 
the quasiparticle mass increase in specific heat measure- 
ments below 10 mK, we have reanalyzed previous magne- 
tization results [23] to estimate C cr in a similar way us- 
ing the Maxwell relation by Paulsen et al. and Sakakibara 
et al. 
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_ 136(] . The static magnetization has shown satu- 
ration below a peak temperature Tp lf23ll . and the normal- 
ized magnetization can be scaled by the normalized temper- 
ature. On the basis of the temperature dependence of mag- 



netization, as given by M 



/3{B)(T/T P ) 2 , below 



T < 0.6Tp, a negative (3(B) is obtained at all measured 
fields and (3(B) asymptotically approaches zero with increas- 
ing field. The field dependence of C a is calculated from the 
equation d{C cr /T)/dB = d 2 M/dT 2 = 2/3 /T P 2 . Thus, the 
C CI /T enhancement that corresponds to a decrease in field 
from 6.21 to 0.20 mT is found to be 26±14 mJ/K 2 mol. The 
anisotropic volume change has been reported in the LFL state 
around IK 02611 . However, the temperature region and the 
physical phenomenon we have discussed are different from 
the experiment. Thus, we assume = 3a a - ignoring the 
anisotropy at ultralow temperatures. The extrapolated value 
of a a /T at zero temperature is —8 x 10~ 6 Kr 2 in 0.016 mT. 
As a result, T cr turns out to be —5 x 10 3 . This growth in the 
critical contribution of the Griineisen parameter is evidence of 
quantum critical fluctuation. 

In summary, the ground state of CeRu2Si2 was studied 
by linear thermal expansion and magnetostriction measure- 
ments. We observed a deviation from the LFL state, indicat- 
ing that the spin fluctuation is dominant at ultralow tempera- 
tures. The negative critical contribution of the linear thermal 
expansion and magnetostriction is clear evidence of the exis- 
tence of an additional QCR The dominant spin fluctuation at 
ultralow temperatures is suppressed by the magnetic field, and 
a crossover occurs from the NFL state to the LFL state. 
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